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Abstract 

We investigate the suppression of 7r°-spectrum in a wide range of p± up 60 GeV/c which is 
caused by the energy loss of the gluon and quark jets in quark-gluon plasma. The physical charac- 
teristics of initial and mixed phases were found in the effective quasiparticle model by analogy with 
previous work 5]. The PHENIX data up 10 GeV/c within the limits of precision are described by 
quasiparticle model with decrease of the thermal gluon mass and effective coupling in the region 
of phase transitin (at T — * T c from above). We also take into account the intrinsic transverse 
momentum k± of partons. 

The suppression factor Raa(j>±) shows the weak rise with increase of p± above 4 GeV/c, then 
it reaches smooth maximum at p± ~ 20 GeV/c and then decrease at p± ~ 60 GeV/c again to 
value Raa at p± ~ 4 GeV/c. The factor Raa in this range of p± is changihg weakly if intrinsic 
momentum k± is taken into account. 
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I. INTRODUCTION 



Energy loss of high energy gluon and quark jets in relativistic A + A collisions leads to jet 
quenching and to suppression of hadron spectra and thus probes the quark-gluon plasma. 
Recently the considerable suppression of tt° spectra was found in central Au + Au collisions 
at v/i=200 GeV [l|. 

The energy loss was investigated, for example, in Ref. Q] in various orders in opacity L/X 
(where L is nuclear radius and A is the gluon mean free path). It was shown that opacity 
expansion is strongly dominated by the first term. It was shown also in Ref |3j that the 
finite kinematic boundaries decrease the energy loss as compared to the asymptotic limits. 
Recently the suppression of 7r° spectra in central Au + Au collisions at y / i=130 GeV was 
investigated in Ref Q with accounting of dominant first term and of the finite kinematic 
boundaries. In this work the physical characteristics of initial and mixed plasma phase 
(i.e. the values T , r , r c ) were used which have been found in Ref || on the basis of the 
quasiparticle model and isentropic expansion. We find that the suppression one can describe 
by quasiparticle model with decrease of the thermal gluon mass and effective coupling in a 
region of phase transition (at T — > T c from above). The possibility of hot glue production 
at the first stage was also considered in Ref The gluon density was calculated with 
accounting of plasma expansion and it was shown that plasma is sufficiently thin : the 
average number of jet scatterings n ~ 1.3, so it is possible to use the model of the single 
hard medium-induced scattering. The physical characteristics of plasma at i/i=200 GeV 
was found by analogy with the Ref. ^. We have used the values T c =177 MeV and /i£=29 
MeV for 200 GeV Au + Au central collisions |6]. These values are consistent with data for 
p/p ratio and for the number of net protons (~ 5) at midrapidity 7]. Using the formulas (10- 
18) in Ref. 0],we have found for gluons (jj^j — 1070 and for hot glue the temperature 
T g ~ 373 MeV. We have obtained also the initial gluon density n g (T g ) ~ 24.77 m\. The 

time which is required to achieve the equilibrium for gluons is: 

dN g 0.587 

Tg = dun J n "" — = °' 834 ^ (1) 
The equilibrium temperature T for both quarks and gluons and corresponding time r can 
be found using the formulas (14-23) in Ref. Q- We have at ^/s = 200 GeV: T =224 MeV, 
r — 1.585 /m n ~ 2.26 fm. The corresponding entropy density is so(Tb) = 58.98 ml and total 
entropy is S = S {T )V ~ 6850. We have also at T = T c = 177 MeV: s c {T c ) ~ 16.72 ml and 



To — 4.5/ m n (i.e. ~ 6.4 fm). In Ref. it was shown, that from conservation of entropy and 
number of net nucleons follows, that the massive constituent quarks (m q and m s ) appears 
with decrease of number of degrees of freedom in the presence of octet of pseudogoldstone 
states. In addition with the same effective number of degrees of freedom appears hadrons 
and resonances in hadron part of mixed phase (for SPS and RHIC at yfs = 130 GeV). One 
can show that it is completely fulfilled for RHIC at ^Js = 200 GeV. 

In Sec. |n]we calculate the energy loss AE of the high energy gluon and quark jets with 
energy E in expanding quark-gluon plasma (at Au + Au collisions at \fs = 200 GeV) in 
the dominant first order of opacity expansion (by analogy with the Ref. |4j) We use here the 
characteristics of expanding plasma which are found in effective quasiparticle model (with 
accounting of the hot glue production at first stage). 

In Sec. IIHI we calculate the suppression of tt° spectra in central Au + Au collisions at 
y/s = 200 GeV in wide range of p± (4 < p± < 60 GeV). In this Sec. we do not take into 
account the intrinsic transverse momentum k± ^ of partons. 

In Sec. IIVI we calculate the suppression of ttq spectra taking into account the intrinsic 
transverse momentum of partons. We use here a Monte-Carlo method for calculation of 
invariant cross section of hadron production. We show, that if take into consider the physical 
restrictions, then the parton cross section does not have the divergence at too large k±^, which 
are considered in other works (for example, in Ref. ^]). Therefore there is no necessity to 
introduce a regulator a 2 in the denominators of the parton-parton cross sections. We show 
that this model reproduces the data for 7To spectra at high p± in pp-collisions at y/s = 200 
GeV without introduction of i^-factor ( at < k\ >~ 1.8 GeV 2 ). We show also, that 
suppression factor Raa{p±) is changes weakly in above-mentioned range of p±, if intrinsic 
momentum kj_,i is taken into account. 

In SeclVl- Conclusion. 



II. ENERGY LOSS OF HIGH ENERGY GLUON AND QUARK JETS 

In previous section the physical characteristics of the plasma stage for RHIC at y/s = 200 
GeV were found. The production of hot glue at the first stage is caused by relatively large 
gg cross section in comparison with qg and qq cross section. We use these characteristics 
for investigation of energy loss of the high energy parton jets. The dominant first order 
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radiation intensity distribution for expanding plasma is given in Ref. |3| and correspond 



to formula (J3J) in Ref. 



.a, 



which have the form: 



| max 



J dz p(z) J d 2 k , a s 



dl _ 9C R E 

dx 7T 2 

z |fe|min 

Qmax 

d 2 q _ a l k q Ti „„„ ( k_ ^) 2 ( z - z ®> 







[q 2 + /i 2 (z)] 2 k 2 (k-q^ 2 



cos ■ 



2x(l - x)E 



(2) 



In this formula the value E is the jet energy and Cr is the color factor of jet: Cr is N c 
for gluon and |iV c for quark jet. It is assumed that quark-gluon plasma can be modeled 
by the well separated color-screened Yukawa potentials. In Ref. |4| was shown, that this 
condition is well realized for partons scattering on gluon potential, but this is not realized 
for partons scattering on the quark Yukawa potential. Therefore the estimation of parton 
energy loss AE on the gluon Yukawa potential is the most real in this model. In formula (j2J) 
we have |g| max = a/3 / u(t)_E', where /i 2 (V) = 47ra s T 2 (r). and the kinematic bounds are 
kmax = mm {4-E 2 x 2 , 4 2 x(l — x)} and k 2 ^ = yU 2 (r) for gluons with the light cone momentum 
fraction x. The value z = r is limited in fact by duration of the plasma phase r g < r < t c 
(and r c < Rau)- The value p(r) is gluon density n g at time r along the jet path, which is 
calculated in effective quasiparticle model. The energy loss AE g of gluon jet is determined by 
integration of dl /dx in (J2J) over x. The corresponding integral Iq(E, t) (look a formula ()22|) in 
Ref. 0]) is calculated in above-mentioned finite kinematic bounds by Monte Carlo method 
for values r(T) from r g to r c and for of E up to 60 — 70 GeV. The energy loss AE g is 
determined by formula: 

&E g (E) = ^ / -4^I (E, r)n 9 (r). (3) 

The values Io{t) and %(t) for expanding plasma at = 200 GeV and the values of 
energy loss for gluon jet we show in Table H] for example at E — 10 GeV. The complete 
energy loss AE is calculated by numeral integration over r. We have AE = 3.6 GeV at 
E = 10 GeV, i.e. not great increase in comparison with = 130 GeV. The complete 
energy loss AE g (E) of gluon jet is shown in Figure Q and the relative energy loss AE g /E 
in Figure El For energy loss of quark jet we have AE q E = | AE g (E). That is important 
non-Abelian feature of parton energy loss j^. 
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Figure 1: The energy loss of gluon jet in quark- gluon plasma at RHIC energy (y/s = 200 GeV) 
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Figure 2: The relative energy loss of gluon jet in quark-gluon plasma. 
It should be noted, that these energy loss correspond to the model of phase transition 



with decrease o 
above. In Ref. 



effective coupling strength G(T) and gluon mass m g (T) for T — > T c from 
we already noted, that model with increase of G(T) for T — > T c agrees 
with new lattice data and provides description of baryon and meson spectra, but it give too 
large suppression of hadrons with large p± which disagrees with data. 



III. SUPPRESSION OF PIONS WITH LARGE P± AT y/S = 200 GEV/C WITH- 
OUT ACCOUNTING OF INTRINSIC MOMENTA K± OF PARTONS. 

The energy loss reduces the jet energy before fragmentation, where the jet transverse 
momentum is shifted on the value AE(E). If this effect is taken into account, we should 
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Table I: The physical values for expanding plasma at y/s = 200 GeV and energy loss A.E g at 
E = 10 GeV 



T (MeV) 


r (m n 1 ) 




/j, 2 (m^) 


n„ {m%) 


h 


AE (GeV) 


373 


0.587 


0.269 


24.34 


24.77 


0.097 


2.64 


350 


0.64 


0.268 


21.35 


20.12 


0.095 


2.61 


325 


0.713 


0.266 


18.27 


15.75 


0.092 


2.36 


300 


0.808 


0.263 


15.39 


12.035 


0.089 


1.905 


275 


0.935 


0.259 


12.74 


8.92 


0.084 


1.62 


250 


1.116 


0.254 


10.32 


6.34 


0.081 


1.34 


224 


1.585 


0.246 


8.024 


3.87 


0.08 


1.065 


210 


1.935 


0.239 


6.85 


2.99 


0.078 


0.96 


200 


2.293 


0.231 


6.01 


2.428 


0.076 


0.84 


190 


2.838 


0.22 


5.16 


1.886 


0.07 


0.7 


185 


3.267 


0.21 


4.67 


1.612 


0.063 


0.59 


180 


3.98 


0.19 


4.025 


1.32 


0.05 


0.448 


177 


4.5 


0.057 


1.16 


1.39 


0.00063 


0.0204 



replace the vacuum fragmentation function by the effective one ^-Dh/ C (z*, Q* ) [10| where 

z * = Z J. (A) 

c _ AE(E) ' V*' 

1 E 

The invariant cross section of hadron production in central A + A collisions is given by 



bmax 



E h —±—= J d 2 bd 2 rt A (r)t A (\b-r\)^2 J dx a dx b d 2 k ± , a d 2 k ± , b 

q abed 



X 



x 9A(k±, a ,Q ,r)g A (k±,b,Q , |b - r\)f a/A (x a , Q } r)f b/A (x b ,Q ,|b-r|)x 
x daz* D h/c (z*,Q 2 ) 
di z c 7iz c 

Here t A is the nuclear thickness function, the fcj_ iQ and k± )b are the intrinsic transverse 
momenta of partons, the values f a / A and f b / A are the partons structure function. It is 
usually assumed that distribution g A (k±) has a Gaussian form. It should be noted, that 
intrinsic k± are more important for final hadron spectra for SPS energies. With the increase 

n 

of energy the spectra become flatter and intrinsic momenta are less important [8|. At higher 
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energy the intrinsic momenta k± are included in order to account for phenomenologically 



the next -to- leading order correction 11 1 (instead of K-factor), The accounting of k± change 
weakly the factor Raa- We will take into account the intrinsic momenta of partons in the 
next Section. 

The upper limit for the impact parameter is 6 max ~ 0.632 Ra u for 10% central Au + Au 
collisions. Since the energy of collisions is very high, we will consider for simplicity the 
collisions of flat disks. In that case the integral of overlapping is Taa(P) = -4n- — 262 fm~ 2 . 
The parton shadowing factor Sa for flat disks also should be taken into account. We used 
the shadowing functions Sa(x,Q 2 ) from EKS98 parameterization These shadowing 

functions affected weakly on suppression factor Raa, especially at large p± (at least at 
p± < 20 GeV/c).It was shown also in Ref. jl7|. At more of p± it is possible the additional 
decrease of Raa due to EMC modification of nuclear structure function. We use the parton 
distribution from Ref. j3] and fragmentation function D g and D g from Ref. ^| both in LO 



parametrization. 

We will use the equation (jHJ) in anot 



rer aspect, in order to compare the energy loss 



. 3 and Q. 



with other works, for example with Ref. 1151 and 161. Let us consider at first the hadrons 
production at parton collisions in the absence of the medium. A parton produced initially 
with transverse momentum p± + u fragmentire into hadron with momentum p±. We have 
z = ^_ u , and we use the factorization scale Q = |j, that is Q(u) = p± ^ u . We have here 
Mmin = and M max = ^j- — p±- The invariant cross section for n° production in pp collisions 
for example gluon jet is (without accounting of partons intrinsic momentum): 

d 6 p 2 J P±[p± + up 



o 
i 

k ( u } 



J d£f gg (£)x 1 G g (x 1 ,Q)x 2 Gg(x 2 ,Q), (6) 



l-fc(u) 



where k{u) = PJ ^ U - We introduce here the new variable x\ = ^ P± ^q^ , x 2 = ^^""^ ■ We 
describe the contribution of the different elementary sections in variable £: 

/ = 1 + <t^Z ( 7) 
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./9<ASj g I ^ 2 I - ^3 ^ 4 IT l) ^2 £4 • V ' 



For corresponding contributions into quark-quark cross section for identical f ld and for 
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different quarks we have 

'm - 9 V £ 4 (£- 1) ' e £ 4 7 27£ 2 ' 



f" = \ e + 1 + e - 1(e ~ 1)3 ) - Al (9 ) 

•'99 nVf/lCt 1 N « t4 / T7 «' V I 



fdif = fd + ^ / 1Q N 

Jqq Jqq ' 27^2 V ' 

In the presence of medium the parton losses the additional energy AE and we have from (jlj): 

z*(u, AE) = P± - (11) 

v ' ; p ± + u - AE v ' 



and also: 



«•(») - § - (12) 



The value z^ ax for gluons (i.e. and uj^) in the presence of medium we find from relation 
(01 4iH 9 = e—ae 9 = ^ f° r ever Y value of p± (Ref. jj|). The value z^ ax for quarks is fined 
by substitution of (AE) q instead of (AE) 9 . The value E min which correspond to gluon jet 
for given p± can be found from Figure d (in large scale). In the range E m - m < E < ^ 
we construct the dependence (^-) 9 (z) using Figure 121 For analogous function for quarks 
(^r-) q (z) we use the relation (AE) q = ~(AE) g . From these relations we find the values 
(AE) gA {u). For example, we find for gluons at p± = 8 GeV (with (3 — 5)% precise): 

/ AE^ 9 

\e~ 

The corresponding equation for quarks is (with accounting that 2^ ax ^ z^. 
( AE\ q 4 



0.42z + 0.05 + 0.35(2 - 0.08)(0.665 - z). (13) 



V E 



V 4 

-j = -(0.4^ + 0.049 + 0.3(2 -0.08)(0.835-z)). (14) 



At increase of p± the dependence (AE)/E become linear. For example, at p± = 20 GeV we 
have: 

' AE^ 9 

( AE\ q . 4 ( AE S 



V E 



0.227z + 0.036 (with 0.2 < z < 0.786), 
V 4 / AE\ 9 

J (^) = -(-_j (with 0.2 < 2 < 0.893). (15) 



We have here in variable u: 



(A£)» = 5.26 + 0.036u (with 5.4 < u < 80 (GeV)). (16) 
The fragmentation scale is now: 

p ± + u- AE°(u) 
QgW = 



2 

P± + u-AE<>(u) 

QqW = ^ • 

Let us write down the formula of type © in AA collisions for 7r° production in gluon jet 
with accounting of energy loss in medium: 

E da^_ = 9 y a s 2 (Q(n)) D g (z*(u),Q g (uj) 



77 d 3 p 2 Tri?^ ./ p±(p±+w) 4 l - Agg W 

"min(Px) P±+« 



J d£f gg (£)x 1 Gg(x 1 ,Q(u))x 2 Gg(x2,Q(u)). (18) 



l-fc(u) 

For calculation of complete energy loss we must take into account also quark-gluon and 
quark-quark collisions with corresponding structure function and elementary cross section 
(171011)) and the energy losses (AE) 9,g (u). We take into account also the gluon collisions with 
sea quarks. 

The suppression factor Raa{p±) is the ratio of the sum of invariant cross section of type 
(|18|) to the sum cross section of binary NN collisions of type© without accounting of nuclear 
effects - jet quenching and shadowing. We have calculated the factor Raa(p±) in wide range 
of p± from 4 up 60 GeV/c. In order to avoid too small values at calculation for great p±, 
we multiply the numerator and denominator in this ratio on value p\_. 

It should be noted, that energy loss in medium can not be expressed by simple manner 
through vacuum cross section da ^ p ^ AE ^ (as for example in Q], Q]). The situation have 
more complex character. 

We show in Table HI the results for suppression factor Raa{p±) for flat disks with ac- 
counting of nuclear shadowing function Sa(x,Q 2 ) from EKS98 parameterization We 
take into account in this Table also the intrinsic transverse momentum of partons (look at 
formula in II V|) . The contribution of nuclear shadowing effects in Raa(p±) is small (less 
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Table II: The suppression Rau,Au(p±) of neutral pions at yfs = 200 GeV. The u m m and z max are 
shown only for gluon jet. 



P-L 




"mm 


z 9 


Raa 


4 


3 


0.572 


0.22 


6 


3.5 


0.632 


0.27 


8 


4 


0.667 


0.34 


10 


4.3 


0.7 


0.38 


14 


4.8 


0.745 


0.42 


20 


5.4 


0.786 


0.44 


30 


6 


0.833 


0.42 


40 


6.2 


0.866 


0.38 


50 


6.5 


0.885 


0.32 


60 


6.8 


0.898 


0.25 


70 


7.2 


0.907 


0.22 



than 10%). It was noted also in the work At large p± > 20 GeV this contribution 

increase up to ~ 20% and that gives decrease of Raa due to EMC effect. 

It should be noted, that at large p± the value Raa{p±) is determined quantitatively in 
the main by quark jets — the contribution of gluons into numerator of ratio Raa become 
small (~ 10% already at p± = 10 GeV) 

The dependence Raa{p±) is shown in Figure EH This function have the weak rise above 
p± = 4 Gev/c and broad maximum in the range p± ~ 15 — 30 Gev/c. The PHENIX data 
there are only up 10 GeV/c It is possible to describe these data by effective quasiparticle 
model with decrease of effective coupling at T —>■ T c from above (look at Table U} without 
additional free parameters. The following decrease of Raa at p± > 20 GeV is caused by 
some factors — in particular, by increase of u m i n with increase of p±. 

It is interesting to note, that dependence Raa(p±) in Figure El disagrees with monotonous 
increase at sufficiently large of p±, which was considered for example in the paper [2^. This 
question have interest and for LHC, where also the continuous rise of suppression Raa with 
p± was found We will investigate this question in the fortcoming work. But of course 
the calculation for large p± may be questionable in very near region to boundary of phase 
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Figure 3: a: The suppression factor Raa{p±) for 7r° in central Au + Au at \fs = 200 GeV in the 
range of p± to 15 GeV/c and PHENIX data, b: The suppression factor Raa which was calculated 
in the range of p± to 60 GeV/c. 

space (p± ~ -y/i/2), where overall energy and momentum conservation may be important. 

IV. THE CALCULATIONS WITH ACCOUNTING OF INTRINSIC TRANS- 
VERSE MOMENTA OF PARTONS 

The invariant cross section of hadron production in central A + A collisions with account- 
ing of intrinsic transverse momenta is given by formula (jSJ. We already have mentioned, 
that we will consider for simplicity the collisions of flat disks. In this case the integral of 
overlapping for Au + Au collisions is TaaO 3 ) — -4s~ — 262 fm~ 2 . We describe the Mandel- 
stam variables for elementary parton-parton scattering in term of light-cone variables x±, 



11 



x-2 



k 2 k 2 

s = xix 2 y/s H — ±LJA. _ 2k_ L1 k_ L2 cos(7i - 72), 

XiX 2 \/S 

t = ( Xiy/s e~ y H ^= - 2k ± i cos(7i" 

z c \ XiJs 



u= -— (x 2 ^~s e~ y + -% - 2k ±2 cos( 72 )l • (19) 

Z c \ X 2 y / S ~ J 

Here p± and y are the transverse momentum and rapidity of produced particles, 
cos(7i) = k-LiP- 1 - ; cos(7 2 ) = fc^p^ ■ It is usually assumed that initial k± distribution 
gN{k±, Q 2 ) has a Gaussian form: 



e 7 ^ 



g N (k ± ,Q 2 )= . (20) 

In equation (jSJ) one should restrict the initial transverse momentum k±i < Xi^/s and 
k± 2 < x 2 \fs such that longitudinal momenta of partons have the same signs as their parent 
hadrons. We introduce the new variables t\ andt2^ k±i = tiXiy/s, k± 2 = t 2 x 2 ^fs, where 
< ti < t max and t max < 1. The variables ()19|) one can describe in the following form: 

§ = sx 1 x 2 V(ti,t 2 , 71,72), 

7T 

a S3C I j v 

t = -—2101^1,71), 

« = -T^xzfcfe, 72), (21) 

where = -^t. 

Here we use designation: 

^1,^,71,72) = 1+^2-2^200^(7! - 72 ), 

Mtuli) = l+t{-2t 1 cos{ 11 )i 

Mh, 72) = 1 + 1\ - 2t 2 cos(j 2 ). (22) 
Now for example instead of f gg (£) © we have: 

f (f f f ^ „ \ = 3(^-^2) _ hMn-fo) 2 , 02 - 2 ) 2 , . 

J 33 U,H,^,7l,72j ^ 2 ^4 ^603 + e(^-02)*0? 



We have analogous modification and for other elementary cross section (jSlHUjl. 

2z c 



In the structure function G g (xi)G g (x 2 ) we have the modification: x\ = and 



x 2 - - 
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The invariant cross section for hadron production (for example in gluon jet) in central 
Au + Au collisions with accounting of partons initial momenta k± have the form: 

d 3 P vR\ J p4(l_^( z )) J 11 J 

x\ X 

2z 

2tt 2tt *X 

J d-fi J dj 2 J d£g 1 {t 1 ,£,z)g 2 {ti,t 2 ,j 1 ,7 2 ,£,z) x 

20 2 (t 2 , 72 )z 

2z¥(t 1 ,t 2 ,7l,72)-'Mtl.7lKj_ 

x /ss(6> *i, *2, 7i> 72)^i^2 s2 ^( if: i^2, 7i72)^iG s (si, Q 2 )x 2 G 9 (x 2 , Q 2 )- (24) 
The value zjL^pi) can be found from condition z^,.„ = — , g x , — r, where the values 

max\r±) max p±+u% nin (p±) ' 

u 9 min {p±) (and also ^ aa , which can be found independently) there are in Table HU 

In__this work we neglect of the transverse momentum smearing from the jet fragmenta- 
tion 

The Gaussians g\ and g 2 in (|24|) have the form: 



exp 

9i(ti,£,z) = 



<kj> 



7l<k 2 ± > 



cxp 



<fci > 



02 *2, 7i, 72, £, z) = 2 . (25) 

With accounting of nuclear shadowing the parton distribution G g (x,Q 2 ) must be modi- 
fied: G g (x,Q 2 )S A (x,Q 2 ). 

In order to find the restriction t max (z), one should take into account, that lower limit 
of integral on £ must be > at any values of t and 7. That correspond to condition: 
^minZ > (j>imax x \- But the value \& have minimum at 71 = 72, and value <j>i have maximum 
at 71 = 7T, therefore we have the condition : 2z(l + t) 2 (l — t) 2 > x"(l + t) 2 , i.e. t < 1— ■^/^- 
At mentioned values of angles 7 the lower limit of £ is: £(2, t) = 2 z(i-t) 2 -x* ' anc ^ we nave 
* < ^max, where t max (z) correspond to £(z,* maa ;) = £ max = ^. For set of values of p± we 
calculate by numeral way the function t max (z). 
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For example, for p± = 8 GeV/c we find: 



, , . /0.04 



- 0.2929 exp (y/z - 0.08)(0.0618 - 10.9786,2 + 27.8446,2 2 - 29.0236z 3 + 10.78z 4 ) (26) 

and for p± = 20 GeV/c we find correspondingly 

[OA 

- 0.2929 exp (y/z - 0.2)(0.6808 - 8.0142^ + 15.83562 2 - 14. 1458/ + 4.744z 4 ). (27) 

At calculation of the invariant cross section by formula (|2*3J) we use Monte-Carlo method. 
We use here the factorization scale Q = |j- and fragmentation scale Q* = |pr- In absence 
of nuclear effects — jet quenching and shadowing the equation f!24|) describes the collection 
of binary NN collisions. We have in this case AE = and ^ aa . = 1. For binary collisions 
we use the value < k\ > pp ~ 1.8 GeV 2 Q-Q- 

In the presence of medium we use the estimation in Ref. jlS^ : 

< k\ >aa=< k\ > pp +C(u m - 1), (28) 

where C ~ 0.4, v m ~ 4, i.e. we have the estimation < k\ >aa— 3. At calculation we use 
also the estimation < k\ >aa— 4. This have a weak influence on the value of suppression 
Raa(p±)- 

We use now the formula (|24|) for calculation of the spectrum of high p± neutral pions in pp 
collisions at \fs = 200 GeV. We take into account the sum of elementary parton - parton cross 
section (jTlllO)) with accounting their modification of type (|23|) . We have not the divergence 
at calculation and therefore we do not introduce a regulator /i 2 in the denominators of 
elementary cross section. We do not introduce also the K-factor. In Figure |U we show our 

n 

calculation for the spectrum of high p± neutral pions in pp collisions and PHENIX data [20J 
up to p± ~ 14 GeV/c One can see that calculation by formula of type (|2*lj) for pp collisions 
reproduces the data in the region p± > 4 GeV/c. 

It should be noted, that we calculated also the factor R A A (p±) at \fs = 200 GeV by 
formula (|2*4^) with Cronin effect and shadowing, but without the energy loss (for flat discs) 
using the value < k\ >aa in the presence of medium (J2Zj). We find the values Raa = 1-56 at 
p± = 4 GeV/c and Raa =1.1 at p± = 8 GeV/c. That agrees with results of paper [isl | . We 
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Figure 4: Neutral pion spectra in pp collisions at \/s = 200 GeV. The data from PHENIX and 
the solid line is calculated result. 

use also the formula f° r calculation of suppression Raa for great p± with accounting of 
intrinsic momentum k±j. They give the small increase of factor Raa - less than 10%. The 
results are shown in Table ITT1 and in Figure El 

It should be noted, that formula of type (}2"4"j) can De used also for investigation of the 
Cronin enhancement at SPS energy (-^/s = 17 GeV). 



V. CONCLUSION 



In this work we use the physical characteristics of initial and mixed plasma phases which 
were found in the effective quasiparticle model for investigation of suppression of neutral 
pion spectra in central Au + Au collisions at yfs = 200 GeV. The specific feature of this 
model consist in decrease of thermal gluon mass and effective coupling in a region of phase 
transition. We use also the hypothesis of hot glue production at the first stage of plasma 
expansion. We describe the following evolution of plasma in isentropic model. The energy 
loss of the high energy gluon and quark jets in plasma leads to suppression of hadron 
spectra. For calculation of energy loss we use opacity expansion 3], which is dominated by 
first term. At calculation we take into account the finite kinematic boundaries, which often 
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do not accounted for. We calculate the gluon density in plasma and show that plasma is 
sufficiently thin. This to some extent justify the model of the single hard medium-induced 
scattering. 

We show that suppression factor Raa(p±) has no continuous rise with increase of p± - 
one has the weak rise, then smooth maximum at P± ~ 20 GeV/c and then it decrease with 
further increase at,!. This qU e Stl0 „ is interesting for LHC, where the calculate QQ 



show continuous rise. We will consider this question in the fortcomin work. It is interesting 
to investigate also the hard photon suppression - apparently it is much smaller than for 
hadrons HQ- 

We use in this work the simplified model of collisions of flat discs at y/s = 200 GeV. We 
are planning to use the more complex model of nuclear geometry in the forthcoming work. 
It should be noted, that model of collisions with transverse expansion and Woods-Saxon 
nuclear distribution with some parametrization of energy loss was used in the work 
There was found also slight increase of nuclear modification factor Raa with increase of 
p± at 200 GeV (at accounting of both gluon and quark jets). In recent work we used the 
parameters which were found in effective quasiparticle model and nuclear modification of 
structure function from [l2J . We take into account also the intrinsic transverse momentum of 
partons k± ti . We show, that one can avoid of divergence at large k^^, if we take into account 
the physical restrictions on initial transverse momentum. Therefore there is no necessity to 
introduce a regulator /i 2 in parton-parton cross-section. It can have the meaning at more 
low SPS energies, for example at investigation of Cronin effect, where the suppression factor 
Raa can be much larger than 1. 
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